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Abstract 

We present new supergravity backgrounds generated by N^. D5-branes, wrapping the S"^ of 
the resolved conifold, in the presence of Nf = 2Nc smeared flavor D5-branes. The smearing 
allows us to take their backreaction on the geometry into account. We discuss the consistency, 
stability, and supersymmetry of these types of setups. We find near horizon geometries that 
we expect to be supergravity duals of SQCD-like theories with Nf = 2Nc. From these 
backgrounds we numerically extract rectangular Wilson loops and beta functions of the dual 
field theory for the regime where our approximations are valid. 
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1. Introduction 



The gauge/gravity correspondence [lj,[2j (for a review and additional references, see [3j) is 
undoubtedly a powerful tool to study strongly coupled gauge theories. Within this frame- 
work, the Chamseddine-Volkov-Maldacena-Nufiez (CVMN) [1],[5],[6] and Klebanov-Strassler 
[7] backgrounds are landmarks in our understanding of = 1 Super- Yang-Mills gravity du- 
als. These backgrounds have been thoroughly studied and extended in many ways in an 
attempt to obtain more realistic models. In particular, since they do not contain fundamen- 
tal matter, a natural question that arises is how to include dynamical quarks. The flavor 
degrees of freedom imply the addition of an open string sector. This is achieved by adding 
flavor branes. In the probe approximation [8] a small number of flavor branes is added to 
the background and a new, tunable, scale rrig appears in the problem leading to interesting 
effects in the predicted meson spectra, the understanding of chiral symmetry breaking and 
phase transition [9]. But in order to neglect the backreaction of the flavor branes we are 
bound to take Nf/Nc 0, in the large Nc limit, and thus, some of the SQCD physics will 
not be captured by this approximation. Recently, Casero, Nunez and Paredes [10] proposed 
a procedure to add a large number of flavors, Nf/Nc ~ 1, to the CVMN background. The 
idea is to introduce Nf flavor branes that extend along the same 3 + 1 gauge theory dimen- 
sions as the Nc color branes that generated the background, wrap an appropriate cycle and 
are smeared along the other angular directions^ In [10] the authors proposed this setup, 
solved the BPS equations asymptotically (for small and large values of the radial variable) 
and studied several gauge theory aspects of the solutions. This flrst step in the quest of a 
SQCD dual has opened a plethora of issues to be studied, tested and understood. 

In the present work we flnd new families of supergravity solutions for the particular case 
of Nf = 2Nc. We consider flavor branes extending along all non-compact directions and 
wrapping a trivial cycle in the compact space. We remark on the consistency, study the 
stability and supersymmetry of this embedding, and present several families of solutions. 
We flnd two two-parameter families of solutions and are able to write analytic expressions 
for them in certain regions. We also flnd a one-parameter family of solutions that we can 
only study numerically. The Nf = 2Nc solution presented in Appendix D of [lU] is included 
in our solutions. Our study of the Nf = 2Nc solutions is based on solving BPS equations. As 
such, all solutions presented are supergravity solutions but not all of them are dual to gauge 
theories. The near-horizon geometries that are dual to SQCD-like theories are contained 
within our solutions. We identify them as the above-mentioned one-parameter family and 
calculate the Wilson loops and beta functions for these backgrounds. 

The paper is organized as follows: In Section [2] we briefly review the CVMN background, 
a supergravity background dual to A/" = 1 SYM in four dimensions. We then summarize the 
procedure proposed in [10] to include a large number of smeared flavor D5-branes into this 
background taking into account their backreaction on the geometry. In Section [3] we examine 
the validity of this approach and discuss some subtleties. We then describe in detail a new set 
of solutions to the BPS equations for Nf = 2Nc. In SectionlHwe study fleld theory quantities; 
we calculate the Wilson loops and beta functions of these solutions. We conclude in Section 

^The same technique has been used to embed a large number of flavor branes into the Klebanov-Witten [llj 
and Klebanov-Strassler [12] backgrounds. 
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ini pointing out some issues and extensions that deserve further study. In Appendix [X] we 
analyze the stabihty of a probe D5-brane in the CVMN and in our backreacted background. 
In Appendix [B] we derive the BPS equations for our particular setup following [T3] and show 
that a single probe brane in this backreacted background is kappa symmetric. In Appendix O 
we present the derivation of the analytic solutions. 

2. Review of Previous Work 

In this section we review the CVMN background as well as how to embed a large number 
of flavor branes into this background following [lOj. The CVMN background was found by 
Chamseddine and Volkov [4j,[5j and later interpreted by Maldacena and Nunez [6J as the 
near horizon geometry of a large number of NS5- or D5-branes wrapping an S'^. Modulo 
the usual subtleties involving Kaluza-Klein modes it is dual to pure A/" = 1 SYM in four 
dimensions 0. The construction in [10] allows the embedding of a large number of flavor 
branes into this background. This leads to a new supergravity background conjectured to 
be dual to A/" = 1 SQCD with a large number of flavors. 

2.1. The Chamseddine- Volkov-Maldacena-Nunez (CVMN) Background 

In two nice papers Chamseddine and Volkov found a monopole solution to SU (2) x SU (2) 
gauged supergravity in four dimensions [1] and lifted it to a solution of ten-dimensional 
supergravity [5]. The significance of this solution was not understood until Maldacena and 
Nunez rederived it as the near horizon geometry of NS5-branes whose worldvolume was com- 
pactified on an S"^ and noticed that as such it should be dual to pure A/" = 1 SYM in four 
dimensions at low energies [6j|^ 

This can be understood as follows. As is well known, in the UV, the correct description of 
parallel NS5-branes is given by little string theory (LST) [l6],[T7j (for a brief review and 
additional references, see [IB]). In the IR, the S-dual picture involving parallel D5-branes 
is more appropriate [19]. Their low energy degrees of freedom are described by (5 + 1)- 
dimensional SYM with sixteen supercharges. If two of the directions of the brane wrap a 
compact manifold, then, at low energies, the degrees of freedom of the D5-branes are de- 
scribed by a (3 + 1) -dimensional gauge theory. However, if this 2-manifold is curved, in 
general supersymmetry is not preserved. In the CVMN solution, the 2-manifold is chosen to 
be a spher^. So in order to preserve some supersymmetry, the normal bundle of the brane 
has to be appropriately twisted [5T]. For the twist chosen, four supercharges are preserved 
and the low energy degrees of freedom of the branes are described by pure Af = 1 SYM in 
four dimensionsn 



Notice again that this supergravity solution will only be dual to a four- dimensional field 
^This was further supported by the more detailed analysis in [T4].[T5] 

■^Other topologies for the 2-nianifolds are possible and are dual to A/" = 1 SYM with matter in the adjoint 
representation PU . 

^In our discussion we have been somewhat careless about the fields describing the motion of the brane 
in the transverse directions. These acquire a mass due to the twist and disappear from the spectrum at low 
energies. 
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theory at low energies {i.e., for small values of the radial coordinate). At higher energies, 
the modes of the gauge theory start to explore the S"^ and the theory first becomes six- 
dimensional A/" = 1 SYM and then, at even higher energies, the blowing-up of the dilaton 
forces us to S-dualize and a little string theory completes the model in the UV. 
The corresponding supergravity solution, after lifting it up to ten dimensions, has the topol- 
ogy of M^'^ X X S^U The metric in the Einstein frame reads. 



l -2 I „2h/jn2 I „-2nj,^2\ , j„2 , ^ /„,,a Aa\2 



■dxls + (de' + sin' edif') + dr' + - A" 



(2.1) 



where (j) is the dilaton. The angles are spherical polar coordinates 6 G [0, vr] and if G [0, 2it) 
parametrizing a two-sphere. The w°' are the SU{2) left-invariant one-forms on the 5''^, 

w'^ = cos iJjdO + sin -ip sin 9d(^, 

w' = — sinip d9 + cos ip sin 9d(p, (2-2) 
w"^ = d^jj + cos 9d(p, 

where 9, (p, ip are Euler angles on the 3-sphere with conventions chosen such that < < vr, 
< if < 2n and < ijj < Air. The twisting is achieved by choosing an embedding of the 
nontrivial U{1) part of the spin connection into the R-symmetry group 5*0(4) ~ SU{2)pi x 
SU{2)l (the group of isometries of the S^)^ The connection one-forms A" (a = 1,2,3) can 
be written in terms of a function a(r) and the angles {9, ip) as follows 

A^ = -a{r)d9, A' = a{r) sin 9 d^, A^ = -cos9d^. (2.3) 

For the metric ansatz (12.11) . one obtains a regular supersymmetric solution when the functions 
a(r), h{r) and the dilaton 0(r) are 

2r 

= —TV' ^2.4a) 
smh 2r 

e^M-) =rcoth2r 5 -, (2.4b) 

sinh 2r 4 

g-2</,(r) ^g-2^0^^^ (2.4c) 

sinh 2r 

where 0o is the value of the dilaton at r = 0. Near the origin r = the function e'^ behaves 
as e'^ ~ and the metric is non-singular. The solution of the type IIB supergravity requires 
a Ramond-Ramond three-form F(3) given by 

■F(3) = -\{w^ - A^) A{w'' - A^) A{w^ - A') + j^F'' A{w'' ~ A''), (2.5) 



^At first sight the metric suggests B}'^ x M x S*^ x S'^ topology. A more careful study shows that the 
singularity at r = is merely a coordinate singularity and that the solution indeed has M}-'^ x R"^ x 5^ 
topology. 

^A different embedding is possible leading to a theory with JV = 2 supersymmetry [20] , 
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where F'^ is the field strength of the SU{2) gauge field A'^, defined as 



pa 



dA 




(2.6) 



The different components of F"- read, 



F 



a'dr A dO, = a' sin Odr A dip, F^ 



(1 - a^) sin 6'd6l A dip 



(2.7) 



where the prime denotes the derivative with respect to r. 

Let us stress that this configuration is non-singular. We should point out that there is a 
one-parameter family of singular non-abelian solutions that interpolates between this regular 
non-abelian solution and the singular abelian solution of pj with a{r) = 0. This family was 
first found in [T5] . 

2.2. Adding Flavor to the CVMN Background 

A possible way to add fiavor to the CVMN background is to consider supersymmetric em- 
beddings of fiavor D5-branes that extend along the Minkowski directions, along the radial 
direction, as well as along a trivial cycle inside the compact directions [13]. The analysis in 
[I3j is done in the probe approximation along the lines of the original work by Karch and 
Katz [8], i.e., the number of fiavors is taken to be much smaller than the number of colors, 
Nf <C Nc, and thus the backreaction of the fiavor D5-branes can be neglected. 
Casero, Nunez and Paredes [TO], went one step further and presented solutions that incor- 
porate the backreaction of the fiavor D5-branes. The geometries they construct depend on 
the ratio x = Nf/N^ which can be kept of order one; they are singular at the origin, but the 
singularity is a "good" one in the sense of the criterion in [22] which means that the metric 
component gu remains finite in the limit r — 0. Everywhere else the geometry is smooth and 
the curvature small as long as QsNc ^ 1 as will be discussed below. These backgrounds are 
conjectured to be dual to A/" = 1 SQCD in four dimensions with a large number of fiavors, up 
to the same caveats concerning the decoupling of the KK modes that were already present 
in the discussion of the original CVMN background without fiavor. 

The general strategy is the following [TOj: One introduces a deformation of the CVMN back- 
ground due to the presence of fiavor D5-branes, derives the corresponding BPS equations 
(see Appendix B of [TO]), and finally attempts to solve them. The fiavor D5-branes are taken 
to extend along the (x°, x^, x^, x^, r) directions and are smeareclll over the {9,(p,9,(p) di- 
rections. These branes can be shown to preserve the same supersymmetry as the background 
for arbitrary values of the angles 6, (p, 6, (p jl3|. Moreover, the smeared fiavor branes will be 
sources for the RR 3-form, resulting in RR fiuxes in the deformed background that can be 
observed as a "violation" of the original Bianchi identity. 

To this end, one introduces an ansatz for the deformation of the CVMN background (-f^(3) = 
^The process of smearing will be explained in detail in section [31 
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0, F(5) = 0) 



dx? 3 + dr' + e'^^''\dd' + sin^ 66^') + 

2k{r) -| 

(2.8) 



2g(r) 2k(r) 

+ ((t(7^ + a{r)def + (u;^ - a(r) sin6'rfv3)^) + ^^(ly^ + cosOdip) 



where we have set a'gg = 1 and A^^^ has been absorbed into e'^^, e^^, e^'^ and dr'^. The dilaton 
is given by = 4/. The left-invariant one-forms w"" on S"^ are the ones given in (12.21) . The 
RR 3-form field strength reads 

F(3) = ^ [-{w^ + h{r)dd) A (w^ - h{r) sinedip) A (w^ + cos 9di^) 

+b'dr A (-de Aw^ + sin ^^(^ A u;^) + (1 - ^(r)^) sin A ci(^ A w^] , (2.9) 
and automatically satisfies the Bianchi identity (i-F(3) = 0. 

It is convenient to introduce a basis of vielbeins for the above metric ansatz (12.81) 
e^" = e^dx'', = le^+9(^w^ + a{r)de), 

eP = dr = ef~^^dp, e' = \ef^^{w'' — a{r) sinOdip), to ^^^'\ 

eO = ef+'^de, = ie^+^' (w;^ + cosedif), ^^'^^^ 

= e^+^smOd^. 

The new radial coordinate p is related to r hj dp = e~''^^^dr and turns out to be convenient 
in subsequent computations. Written in this basis, the RR 3-form field strength becomes 

F(3) = -2Nce-^^-''^-^e^ A e' A + ^fe'e'^^-^-'^-'^e^ A (-e^ A + A e^) 

^ ^g-3/-2/.-fe(^2 _ _^ ^^gS A A (2.11) 



+ N^e~^^'3'^-\b - a)(-e^ A + A e^) A 



To be able to account for the backreaction of the flavor branes, the action has to be augmented 
by the DBI and Wess-Zumino actions for the flavor D5-branes. The complete action then 
reads 

•S* = iSgrav + 'S'flavor) (2.12) 

where, in Einstein frame, we have 

S,r.. = ^ J d^'x^/^) {r - \{d,cl>){d^cl>) - l^e^Ff,^ , (2.13) 

and 

^flavor = T,y^i- ( d'xe't>''^^) + I P[C(6)] ) . (2.14) 

One of the effects of smearing the Nf — > oo flavor branes along the two transverse 2-spheres 
is that there will be no dependence on the angular coordinates {6, (f, 6, (p) of the functions 
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f(p),h{p),g(p) and k(p) that determine our metric ansatz (I2.8|l l^. significantly simplifying 
the computations. After the smearing we can write 

5flavor = ^(-y" Asin^sin^>/2y^+ j Vol(3^4) A qe)) , (2.15) 

with the definition Vol(3^4) = sin 6 sin d6 A dip A d6 A dip. Once the smeared flavor D5-branes 
are incorporated into the background, the Bianchi identity for F(3) (which is identical to the 
e.o.m. for C(6)) gets modified to 

rfF(3) = ^ sin ^ sin OdO A dcp A dO A dip, (2.16) 

as a result of adding a Wess-Zumino term. This can be solved by adding the following term 
to the original F(3) in (12.91) . 

™ = sin ede AdipAw'^ = __te-3/-2'^-'=e^ A A e\ (2.17) 
The full RR 3-form field strength now reads 

F(3) = -2iV,e-=^^-25-^ei A A + ^fo'e'^z-g-h-fc^r ^ ^_^e ^ gi + ;\ 

^ ^g-3/-2/.-fe(^2 _ 2a6 + 1 - x)e' A A (2.18) 



+ N^e-^f'3'^-^{ij _ a)(-e^ A + A e^) A e^ 

where the only modification is the appearance of the term involving x = ^ in the second 
line. 

The first order BPS equations for this ansatz (see Appendix B, [10]) can be partially solved, 
yielding 

h{p) = (2.19a) 
smh(2p) 

^ iVc 2&cosh(2p)-2 + x 

2 acosh(2p) - 1 ' y ■ J 

e^'' = — (2acosh(2p) - 1 - a^), (2.19c) 



^The validity of this will be discussed in Section [31 
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leaving us with (coupled) differential equations for a(p), k{p) and /(p) 



2 / 2e^'' - Nf (acosh2p- 1)^ . , \ , , 

d.a = ; + sinh 2p -2ap] , 2.20a 

2pcoth2p V '^Nc-Nf sinh2p h hj. \ > 

2 ( 2e^^ + jV 

9o/c = ; ; ( — —^a sinh 2p(a cosh 2p — 1) + 

(2pcoth2p)(l -2acosh2p + a2) \2A/^,- AT/ ^ ' 

2 sinh 2p 



dpf 



+2p(a^ cosh 2p - 2a cosh 2p + 1) ) , (2.20b) 

2p 



-1 + a cosh 2p) sinh ^2p 



4(1 + a2 - 2a cosh 2p)(-l + 2pcosh2p) 



X ^-4p + sinh 4p + 4ap cosh 2p - 2a sinh 2p - ^-^ — asinh^2p^ . (2.20c) 

Equations (12.201) can be solved numerically for x 7^ 2. A detailed discussion of the solutions 
along with asymptotic expansions can be found in [in],[23|- Note that the equations (12.201) 
are not valid for x = 2. In this case (cf. Appendix D of flU]) equation (I2.19ap implies that 
in the limit a; ^ 2 we have h{p) = 0. Furthermore, the authors in observe from their 
numerical solutions that a(p) asymptotes to 1/ cosh2p in the limit x — > 2. In the Appendix 
[B]we derive the BPS equations for 6(p) = 0. For 



1 



cosh 2p 



f2.211 



they reduce to 

e'^-s = ^ tanh 2p , (2.22a) 

= e^^ - N, , (2.22b) 
dpk = -{e^^ + N^)e-'^-^ + 2coth2p, , (2.22c) 
dpcj) = Ad J = AT.e-^-^ . (2.22d) 

Solving this system of equations will be the topic of Subsection 13.21 below. 
As an aside, let us also briefly discuss yet another solution to the BPS equations fl2.19p and 
f l2.20p for X = 2 that was found in [TO] by embedding a large number of flavor branes in the 
singular CVMN background and will be important later. As discussed above, from equation 
fl2.19ap it follows that b{p) = and it is consistent to set a(p) = as wellj^ In pi)] the 
solution to the BPS equations is found to be 

e'' = j, ^'' = i^' ^'' = ^- = 4/ = 0o + P, (2.23) 



^Thc two classes of solutions could be named abelian for a{p) — and non-abelian for a{p) — l/cosli2p 
in analogy with the abelian and non-abelian CVMN backgrounds. We expect that there should again be 
a one-parameter family of solutions interpolating between the two, analogous to the one-parameter family 
interpolating between the abelian and non-abelian CVMN background found in [TS] . 
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where < ^ < 4. The corresponding Einstein frame metric has a singularity for p 0. 
However, this singularity can be smoothed out by turning on a temperature in the dual field 
theory, thus finding a black hole solution. The black hole temperature turns out [24] to be 
reminiscent of the Hagedorn temperature in little string theory. Therefore, it seems that this 
solution is not dual to a field theory but to a little string theory. 



3. New Supergravity Solutions 

In this section we consider the inclusion of a large number of flavor branes into the CVMN 
background for the special case of Nf = 2Nc. We discuss the validity and shortcomings of 
the smearing approach in Subsection 13.11 before solving the BPS equations fl2.22p in Subsec- 
tion [321 In Subsection 13.31 we discuss properties of our solutions. 

3.1. Including a Large Number of Flavor Branes 

We take the flavor branes to extend along the x^,ilj,p directions and place them at fixed 
values in the 6, ip, 6, (p directions. The reason for taking the branes to extend along the non- 
compact p direction is that we want to send their 't Hooft coupling to zero. For our fiavor 
branes we have 

2(27r)^v/^^ 
- Vol{S^) ' ^^-^^ 

where Voli^S^^ is the dimensionless volume of the ^/^-direction wrapped by the branes. For 
energies small compared to the inverse size of the this coupling determines the interaction 
of the 5-dimensional gauge fields living on the fiavor branes and the quarks living on the 
4-dimensional intersection of the color and fiavor branes0 When taking the near horizon 
limit we take a' ^ to decouple open and closed string modes while keeping gs^c large 
but finite to be able to trust the supergravity approximation. In this limit the interaction 
between the quarks and the five dimensional gauge fields disappears and we obtain SQCD 
with a global fiavor symmetry. Note that this would not be possible if the fiavor branes 
wrapped a compact 2-cycle. 

As mentioned above, we have to include the DBI and Wess-Zumino actions f l2.14p for the 
fiavor branes. Using delta functions these can be written as integrals over the 10-dimensional 
spacetime. However, solving the equations of motions containing delta functions is hard. So 
we invoke a 'smearing' process first proposed in [25j. If this smearing process is to preserve 
supersymmetry, we have to ensure that all the branes at different positions in the transverse 
space preserve the same supersymmetry as the background. This was checked in [13] and can 
easily be seen from the results in Appendix [B] Formally, this smearing process can be thought 
of as follows: As is well known the delta function can be written as an infinite sum using 
the completeness relation for the eigenfunctions of the Laplacian on the space transverse to 
the branes. The smearing corresponds to truncating this sum to the zero mode. Intuitively, 



^°It should be noted that the symmetry is formally a global symmetry due to the presence of a non-compact 
direction but there still is a non-vanishing interaction between the quarks due to their charge under this 
symmetry unless A/ = 0. 
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however, one might want to evenly distribute a very large number Nf of branes over their 
transverse directions, i.e., in our case the 9, ip, 9, if directions. It is not immediately apparent 
that these two prescriptions are equivalent. To check this, one should make an ansatz for all 
the fields in the type IIB action, the DBI action, and the Wess-Zumino action that depends 
on the positions of all the branes. Then one can expand all fields and the delta functions 
in two sets of spherical harmonic functions corresponding to the two sets of coordinates 9, ip 
and 9, ip. The hope then is that there exists a consistent truncation to the zero modes and 
that in the large Nf limit the contribution of the higher modes is negligible. This seems 
plausible due to the homogeneous distribution of a very large number of flavor branes. In 
this paper we do not attempt to verify this hope but rather try to provide some checks that 
indicate that the interpretation of our solutions as being dual to A/" = 1 SQCD might be 
consistent. Note that our solutions are supergravity solutions with Af = 1 supersymmetry 
since we solve the BPS equations. However, to conjecture that they are dual to a gauge 
theory would require that the smearing prescription provides a valid approximation to the 
brane configuration. 

Before we move on, some more comments on the consistency of our setup are in order. The 
flavor branes extend along all non-compact directions. Gauss's law then implies that these 
branes cannot carry a net charge. As in [8], we accomplish this in our setup by wrapping 
the branes around a topologically trivial cycle in the compact space that carries no flux. To 
be more specific, our flavor branes wrap the Hopf fiber, i.e. an inside the S^, which is 
clearly trivial. Most stable D-brane configurations are stable precisely due to the charges 
they carry. So we should ask if these flavor branes will not simply slip off the S^. A stability 
analysis shows that the slipping mode for a probe brane embedded in this way into the 
CVMN background indeed has a tachyonic potential as expected but is in fact stabilized by 
effects similar to the ones observed in [26],[27j. In other words there exists the analogue of 
a Breitenlohner-Freedman bound in this geometry. Furthermore, one can show P3] that the 
embedding is kappa symmetric, which implies that the configuration we would like to smear 
is not only stable but also supersymmetric p8]0 For the interested reader the stability 
analysis is carried out in Appendix \M 

To verify that the new background in which the backreaction of the smeared flavor branes 
is taken into account still preserves the same amount of supersymmetry as the CVMN 
background, we show in Appendices \M and [B] that an additional probe brane extending 
along x^,ip,p and at fixed 9,ip,9,ip is still both stable under small fluctuations and kappa 
symmetric implying that the system is indeed supersymmetric |28j . 

Since we are placing our flavor branes at different positions in the transverse space they 
support a U{l)'^f and not a U{Nf) gauge group. However, at energies small compared to 
the inverse size of the compact transverse space the theory will not be able to resolve the 
transverse space and should effectively behave as a theory with U{Nf) symmetry. So our 
compact directions need to be big enough to trust supergravity but also small enough so 
that, at the energies we are interested in, we can neglect KK modes and have a U{Nf) flavor 
symmetry. 

^^The fact that the brane probe is a solution of the worldvolume equations of motion together with kappa 
symmetry aheady impHes supersymmetry. 
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3.2. Solutions of the BPS Equations for Nj = 2Nc 

In this subsection we present new supergravity solutions describing Nj = 2Nc smeared flavor 
branes in the CVMN background, extending the result in Appendix D of [TO] . 
Our metric ansatz is (12.81) 

ds% = e^f^P^ (^dxl, + e^'^^^^dp^ + e^^^P\de^ + sin^ Od^^) 

+ {{w^ + a{p)def + {w^ - a{p) sin 9d^f) + ^^{w^ + cos ed^Yj , (3.2) 



cosh 2p 

is 



where a(p) = ^^J^ and the w"" are given in (12.21) . The RR 3-form in the vielbein basis (12.101) 

2 

+ N^a e-3/-f-/^-fc(e« a + A e^) A e\ (3.3) 



The BPS equations are derived in Appendix [B] and reduce to e'^ ^ = |tanh2p, 

dpe''+^ = e^^ - N„ (3.4a) 

dpk = -{e^^ + N^)e~^^^ + 2coth2p, (3.4b) 

and 

dpct) = Adpf = N,e-'-\ (3.5) 

Note that the dependence of the solution is trivial since rescaling e^'^ and e^"'"'* by a factor 
of Nc eliminates from the equations. 

Once a solution to the coupled first order equations ( 13. 4p is found the dilaton (j) can be 
obtained by integration of (13.51) . The equations (13. 4p do not seem to allow us to write the 
general solution in a closed analytic form. The equations can easily be solved numerically 
and there are in fact a number of things that can be done analytically. We can find an 
analytic solution valid for small values of p in the form of a power series. To leading order 
the solution takes the form 



e^+VA^, = po-p, (3.6) 

'p_ 
.Pi 



e''/K = (-) , (3.7) 



implying that, at least for all solutions of this form, e^^ goes to zero for small p. 
Instead of studying the small p behavior we can find closed analytic expressions approximat- 
ing the two-parameter family of solutions to the equations (13.41) in the following regions: 

Region I: e^'' < 

Region II: 6^+^^ > N^, e^^ (3.8) 

Region III: e^^ > A^^, 
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1 2 3 4 5 6 Nc 



Figure 1: This plot shows examples of numerical solutions of the above system of equations 
( 13.41) as well as the regions in which the equations can be solved analytically. We plot three 
solutions starting at both e^^^ /Nc = po = 1/2 and e^^^ /N^ = po = 2 and one starting at 
e^+ViVc = Po = 5. 

These regions, together with some examples of numerical solutions, can be seen in Figure [TJ 
There are two qualitatively different two-parameter families of solutions. For small val- 
ues of p both start out at e^*^ = 0, e^'^'^/Nc = po- For large p one of the families crosses 
e^'^ = |e^'^'' — ^ and then approaches it from above. The other family reaches a singularity 
at e^'^ = e^~^^ = for a finite value of p. For po < Pc ~ 1-6, we can make a further distinction 
among the singular solutions. For given po < Pc and large enough values of pi the func- 
tion e^~^^ decreases monotonically while below a certain value of pi the solutions describe 
a loop in the e'^^-e^^'^ plane. A schematic plot showing the different regions in the po-p/ 
parameter space corresponding to the three classes of solutions is shown in Figure [21 The 
singular and non-singular two-parameter families are separated by a one-parameter family 
that approaches e^'^ = e^~^^ = Nc for p ^ oo. In fact, it asymptotically approaches the ^ = 2 
solution found in [TO] and given above in equation (12.231) . This one-parameter family will 
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Figure 2: This plot shows the regions in the po—pi parameter space corresponding to different 
classes of solutions. Both shaded regions correspond to the solutions with a singularity at 
finite p. For the region shaded in lighter gray e^^^ is monotonically decreasing. The region 
shaded in darker gray corresponds to the solutions describing a loop in the e^^-e^^^ plane. 
The white region corresponds to solutions that approach e^'^ = |e^+'* — ^ for large p. 

play a special role as discussed in the next subsection. 

The analytic solutions discussed below comprise two-parameter families originating from Re- 
gion I as long as they remain inside Regions I, II or III for all p as well as a one-parameter 
family contained entirely in Region III. We will merely summarize these analytic solutions 
in this section and refer the interested reader to Appendix O for a detailed derivation. 
Unfortunately, we were not able to find an analytic solution for the one-parameter family 
separating the two two- pa rameter families. We will therefore have to study them numerically 
in the next subsection. 
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The Solution Contained in Region I 

For e^^ -C A^^^ we can solve the BPS equations (13. 4p to first order in e^^. The solution is 



e'+ViVc = PO-P + 



gj^(i(l+8(p-p„fsmh4p) 



- (p - Po) cosh 4p - Po + gPlp^ - 3ppo + 3po) j , (3.9) 
e-/iV, = {^{i-IL]]'. (3.10) 



V 2p7 V Po 



^^One particular solution of the one-parameter family of solutions ending at the point e^*^ = e^+'' — Nc 
and coming from the overlap of Region I and II was already found in Appendix D of [10 . 
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For this solution to be entirely contained in Region I we need to impose e^*^ ^ Nc for all 
p. Note that this solution does not only have a singularity at p = but also a second one 
near p = po- At the singularity at p = the dilaton and the metric component gu both go 
to a constant. On the other hand both go to infinity at the other singularity. It is not clear 
whether these solutions can be given a physical interpretation. 

The Solution Starting in Region II and Ending in Region III 

For a somewhat limited choice of initial conditions, meaning po ^ 1 and p/ ^ 1 or po, p/ ^ 1, 
we can find a closed form for the analytic solution approximating the solution starting in 
the overlap of Region I and II and ending in Region III. For p/ <^ 1 and po ^ 1 we find 

while for ie^^" 3> p/ 3> 1 and po ^ 1 the solution is 

e^+VA^, = Po l + TH ? -P, (3.13) 



64poP/ 
sinh2p\ ^ / 3e^P 



As discussed below, we think that these solutions describe the geometry generated by the 
color and fiavor D5-branes. For somewhat less restrictive initial conditions the matching can 
still be done but with piecewise defined functions that give a slightly better approximation 
to the numerical results. This is described in more detail in Appendix [Cl 

The Solution Contained in Region III 

Solving the equations for e^^ ^ Nc by neglecting Nc in the BPS equations (13.41) one finds 

e^+^/Nc = A(85- 12p + 3sinh4p)^ , (3.15) 

e^'/Nc = 8Asinh^2p 

(85-12p + 3sinh4p)3 

For this solution to be entirely contained in Region III we need to demand that e"^^ ^ Nc 
for all p. In particular this has to hold in the limit p which implies B = 0. After 

4 

making the identification A = 4^1 this solution can be seen to correspond to the deformed 

2333 

conifold [29] 

e^-^ViVc = ^^(^^"h4p-4p)^ ^ 

23 

-J (smh4p — 4pj3 



14 



In order for it to be contained entirely in Region III we require e ^ (|) * ^ 1. 
3.3. Discussions of the Solutions and the Decoupling Limit 

The two-parameter families of solutions reaching the point e^^ = e^~^^ = for finite p = ps 
are singular both at p = as well as at p = p^- The Ricci scalar close to the singularities 
scales as 1/p^ and l/(p — Ps)^, respectively. We do not know whether these solutions are 
physical. 

For the two-parameter family of solutions which approaches e^'^ = |e^+'^ — ^ as p — >■ cxd we 
can write the metric for large p as 

ds' = a'^.iV,e<^W/^ f 4% + + ^'^4nl ^ (3-19) 
\a'gsNc "y 

where we defined a new radial coordinate r ~ e^'^^^ and ds'^_^_^ is the Einstein metric on the 
homogeneous space {SU{2) x S'f/(2))/[/(l)diag- So for large values of the radial coordinate 
the metric approaches a product of Minkowski space and the singular conifold. As p goes 
to zero the geometry is singular and the Ricci scalar diverges like 1/p^. This should be 
compared to the usual 1/p divergence for fiat D5-branes. We attribute the discrepancy to 
the way the flavor branes are embedded. According to the criteria developed in |30j the 
solutions should still be accepted as being physical since the gu component of the metric 
remains bounded as p approaches zero. The dilaton increases monotonically with p and for 
large values of p approaches a constant. So as long as gsNc 3> 1 the Ricci scalar is small for 
all large enough p, and, if (7^ ^ 1 as well, this means that the supergravity approximation 
is valid. We believe that these solutions describe the geometry generated by the color 
branes in the presence of Nf = 2Nc smeared flavor branes. 

These two two-parameter families are separated by a one-parameter family that approaches 
g2fc _ ^g+h _ p qq_ This situation of having two classes of solutions separated by 

'special' solutions should be familiar from e.g. the Klebanov-Tseytlin solution |31j. There 
one has an integration constant Bq. For 60 < one flnds solutions that have a second 
singularity for flnite value of the radial coordinate. For 60 > the solutions describe the 
geometries generated by the brane conflguration. For Bq = one obtains special solutions 
describing the corresponding near horizon geometries. Before flxing the boundary conditions 
at inflnity a similar behavior can be observed for a stack of Dp-branes in flat space. This 
leads us to the conclusion that our one-parameter family separating the singular and non- 
singular two-parameter families describes the near horizon geometry for our brane setup and 
should be dual to A/" = 1 SQCD with Nf = 2Nc. Due to the singularity at p = we should 
impose an IR cut-off. There also has to be a UV cut-off due to the presence of Kaluza- 
Klein modes. These become important as the energy increases. At even higher energies the 
theory becomes six dimensional and has a little string theory as its UV completion. We 
should note that it is not clear that the IR and UV cut-off are even compatible (in the sense 
that all approximations are valid) but we flnd that at least for a certain range of initial 
conditions this is the case. In the next section we study Wilson loops and beta functions 
for these solutions to support our interpretation. Unfortunately, we do not have an analytic 
expression describing these geometries because the near horizon limit takes us out of the 
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regime where our approximations (13.81) are valid. So we study these solutions numerically. 
Figure [3] shows a few examples of this one-parameter family of solutions. 
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Figure 3: This plot shows representatives of our one-parameter family of solutions sepa- 
rating the singular and non-singular two-parameter families starting at e^^^ /N^ = Po = 
0.1, 0.5, 1, 1.4, 2 and e^'^ = 0. Longer dashes correspond to larger po values. 
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4. Field Theory 



In this section, we will use the gauge/gravity duality to study rectangular Wilson loops as 
well as beta functions of the field theories dual to our supergravity solutions. 

4.1. Wilson Loops 

As usual the Wilson loop for two non-dynamical quarks can be computed by evaluating 
the Nambu-Goto action of a fundamental string with appropriate boundary conditions on 
shell [32], [33]. We will parametrize the string world sheet hj t = t, x = a and p = p{cr)- 
Suppose we study a quark and an antiquark separated by a distance L along the x-direction 
in gauge theory coordinates. The energy and length of the string can be given as functions 
of p, where p is the minimal coordinate distance between the worldsheet and the singularity. 
In the string frame these become 

E(p) = ( 2 / dp' - 2 / e*(^')+'=(^')dp' ) , (4.2) 

where p is a cutoff that can be taken to infinity at the end of the calculation. We will study 
these integrals numerically and must work with a large but finite cut-off. We have checked 
that the cut-off is large enough to make our results cut-off independent. The results can 
be used to obtain a relation E{L) = Vgg{L). We will calculate the Wilson loops for the 
same set of parameters as shown in Figure [3] above. For solutions with po < Pc ~ 1-6 we 
find that L{p) has a maximum at p comparable to the value of p for which the supergravity 
approximation breaks down. This is shown in Figures H] and [51 For po > Pc ~ 1-6 L{p) 
approaches y/ a'ggNc'K from belowj^ We can only trust the results of our calculations in a 
regime where the string worldsheet is entirely contained in regions where the supergravity 
approximation is valid and the KK-modes are negligible. The scale at which the KK-modes 
become important was identified in [TD] as the value of p for which L{p) reaches the value 
L{oo) = TTy/a'QsNc. For all solutions this occurs around p ~ 2.5. The value of p above which 
the curvature is small and supergravity can be trusted depends on the particular solution. 
The Ricci scalar for some representatives can be seen in Figure [51 For po < Pc the curvature 
is small for p > 1.2 while for large po we can only trust supergravity for p ^ Po- Note 
that this means that for po > Pc ~ 1.6 there is no regime in which both supergravity can 
be trusted and the KK-modes can be neglected. For small po we can trust our results for 
1-2 ^ P ^ 2.5. The quark- antiquark potential extracted from these Wilson loops for this 
range of the radial variable is shown in Figure [61 Over the range of radial variable p for 
which our approximations are valid the function L{p) is a monotonic function of p and we 
do not observe any signs of pair creation. We should not try to interpret the maximum 
in L{p) in any way since it occurs for values of p where the supergravity approximation is 

^•^Using the large p behavior of our solutions, namely that e^*^ — > Nc and e"^ e'^°'^P, it is straightforward 
to check that for large p we have L{p) ~ , — n^/a'gsNc . 
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Figure 4: This plot shows the separation between the endpoints of the string L as a function 
of the minimal distance between the string worldsheet and the singularity at p = 0. L has 
units \/a'gsNc. We have dashed the lines so that longer dashes correspond to larger po 
values and po takes the values po — 0.1, 0.5, 1, 1.4, 2. We see that L has a maximum for all 
Po < Pc ~ 1-6. 

ff'gsNc R 



1.5 



0.5 



-0.5 




Figure 5: This plot shows a'gsN^R versus p for the near horizon geometries for po = 
0.1,0.5, 1, 1.4,2. Again longer dashes correspond to larger po values. The Ricci scalar has 
a maximum at a p value which is comparable to the value for which L has its maximum if 
Po < Pc^ 1-6. 



no longer valid and one should expect stringy corrections to alter the behavior of L{p) for 
small p. We do think that the dual field theory contains dynamical quarks and their effects 
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Figure 6: The quark-antiquark potential for po = 0.1,0.5,1,1.4. E is given in units of 
^'''2^^" ^^"^ ^ units of y/a'gsNc- 



should certainly be observable on the supergravity side but we think that pair creation in 
our models occurs in a regime where stringy corrections can no longer be neglected. One 
might speculate that L{p) remains monotonic for all values of p once stringy corrections are 
taken into account and that L{p) remains finite as p approaches zero signaling pair creation. 

4.2. Beta Functions 

Recall that the Wilsonian beta function for SQCD is given by 

/?. = -^(3iVc-iV^(l-7o)), (4.3) 

where 70 is the anomalous dimension of the quark superfield. We choose to calculate the 
equivalent quantity 

There are two distinct ways to establish a relation between the Wilsonian cut-off p and the 
radial direction p discussed in the literature [Mj, sometimes referred to as the holographic 
relation and the stretched string relation. Since we only know our solutions numerically, and 
the stretched string relation is much easier to implement, we will restrict ourselves to the 
latter. In other words we use the following relation between the Wilsonian cut-off and the 
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radial distance 

27rv a' Jo 

This is shown for our usual representatives in Figure [3 




12 3 4 



Figure 7: This plot shows the Wilsonian cut-off as function of p for po = 0.1, 0.5, 1, 1.4, 2. /i 
is given in units of ^'^"v^gf^^ and we have dashed the lines so that longer dashes correspond 
to larger po values. 




Figure 8: This plot shows the beta function (14. 4p for po = 0.1,0.5,1,1.4,2. p is given in 
units of '^t^Jvo^'' ^^"^ have dashed the lines so that longer dashes correspond to larger po 
values. 
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The range of 1.2 < p < 2.5 for which both the supergravity approximation is vahd and the 
KK-modes are neghgible corresponds to one order of magnitude in energy for the solution 
starting at e^^^ /Nc = Po = 1-4 and e^'^ = 0. The range can be made arbitrarily large by 
making po small. 

In Figure [H] we show the numerical results for the beta functions. For the range of 1.2 < p < 
2.5 the beta functions are monotonically decreasing. They start at values smaller than 0.5 
and approach zero for large values of /x. This means that the anomalous dimension is always 
negative and approaches —1/2 from above in the UV. 

5. Conclusions and Future Directions 

In this paper we found new supergravity backgrounds incorporating the backreaction of a 
large number Nj = 2Nc of flavor D5-branes. Among the solutions presented we identify a 
one-parameter family of solutions that describes the corresponding near-horizon geometries. 
We therefore expect this family to be dual to an SQCD-like theory. We calculated the 
Wilson loops and beta functions of these backgrounds finding a behavior that is consistent 
with the claim that it is dual to an A/" = 1 SQCD-like theory with Nf = 2Nc matter in 
the fundamental representation. There are several issues related to the present work that 
deserve further study: 

• The field theory dual to the backgrounds presented here needs to be better understood. 
This problem can be approached from different, complimentary, angles. From the su- 
pergravity side it would be interesting to calculate field theory observables like spectra 
of glueballs, mesons and baryons and compare them with results for the CVMN back- 
ground [I3],[35j. Seiberg duahty and its geometric realization needs to be explored. It 
would also be very desirable to have a detailed construction of the field theory from 
the open string side. 

• The universality of the shear viscosity to entropy density coefficient [M| , [2Z] , [3S] was 
a successful result that opened the possibility of studying real life finite temperature 
QCD with gauge/gravity methods. Several QGP signatures have been investigated 
[59] yielding a promising picture. However, all the non-extremal gravity backgrounds 
used to date do not incorporate dynamical quarks. It is unclear if finite temperature 
A/" = 4 SYM can capture the dynamics of QCD above the deconfinement transition 
with enough precision to be useful in RHIC physics. Thus, finding a black hole solution 
of the backgrounds presented here and studying its hydrodynamic properties as well as 
QGP signatures is undoubtedly an interesting problem but we suspect that one might 
find problems similar to the ones encountered in |15j. There it was observed that in 
the CVMN background the SYM modes do not seem to cleanly decouple from the 
massive Kaluza-Klein modes and there seem to be no black hole solutions for which 
the Hawking temperature is less than the Hagedorn temperature of LST. It still seems 
worthwhile to understand whether these problems persist in our setup. 

• Another avenue to explore in the quest for a QCD-like dual is finding a stable non- 
SUSY deformation [10] of the backgrounds presented here. 
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A. Appendix: Stability Analysis 

In this appendix we will discuss the stability of probe D5-branes in both the CVMN back- 
ground and our background that takes the backreaction to the presence of smeared flavor 
branes into account. The probes will extend along the four- dimensional spacetime, along the 
Hopf fiber inside the three-sphere parametrized by ip, and along the radial direction. While 
Euler angles are convenient coordinates on the three-sphere in the rest of our paper, they are 
not very convenient for the stability analysis and we change to spherical polar coordinates 
on the as follows 

v4 e (p + ii) . 

X = cos - cos — ^ — = smi(7cosf, (A.laj 

X^ = cos - sin ^ = cosw, (A. lb) 

2 2' ^ ' 

X = sm - sm = — sm w sm v sm u , (A. Ic j 

. 6 (p — if) 
X = sin - cos — ^ — = sin If sin f cos M . (A. Id) 

If we denote the worldvolume coordinates by (.^^, a, p)0 the cylinder solution we are inter- 
ested in is described by the embedding 

x'^(e,a,p) =e, M(e,cr,p) =a, 
r(^'^,a,p) = p, t;(^^,a,p) = f , 

^(e,a,p)=^o, ti;(e,^x,p) = |, ^ • ^ 

<^(^'',(^,P) = 'ho- 
using reparametrization invariance of the worldvolume theory, we can write the most general 
fluctuations around this background as follows 



x'^(e,cr,p)=e, M(e,a,p) = a, 

r(e, a, p) = p , ^, P) = f + CT, p) , 

^(p, a, p) = ^0 + cr, p) , w{i^ a, p) = f + a, p) , 

^"^We take a to have period 2-k. 
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The fluctuations a, (3 correspond to translations on the two-sphere and it is easy to see 
that their zero-modes will have a flat potential allowing us to drop them from our analysis. 
Furthermore, we can convince ourselves that a and ^* dependence of a and (3 will yield a 
positive contribution to their energy. Therefore, when analyzing the stability, we can limit 
ourselves to fluctuations of the form 

r{,i'',cr,p) = p, 
^(e,or,p) = 0o, 



^;(e,a,p) = | + a(t,p), 
w(e^a,p) = f + /3(t,p). 



(A.4) 



A.l. Stability of the Cylinder Solution in the CVMN Background 

We are now ready to study whether or not there are unstable modes corresponding to the 
brane slipping off the three-sphere. The worldvolume action is given by 

5flavor = -T5/ (f XC*'^ ,/^) + T, I P[Ci^,)l (A.5) 
J Ma J Me 

where (7(6) as usual is the determinant of the pull-back of the metric on the target space and 
P[C(6)] denotes the pull-back of the six-form via the embedding map. The six-form is given 
by 

C(g) = — A dx^ A dx'^ Adx^ AC, (A.6) 
8 

where C is given by 

C = {{a^ - \)o?e-'^^ - IGe^'*) co^Qd^ A dr 

- [a^ - \)e-'^^uf' Adr + a'(sin Odip Au^ + dO A uj'^) . (A.7) 

For our embedding we flnd that, in units of the brane tension, this leads to a potential for 
the fluctuations given by 

V{a,(3) = e^"^ ^cosacos/3- |(1 -a2)e-2^cos2acos2/5^ . (A.8) 

This depends on the radial coordinate via the dilaton as well as the functions a and h. It 
turns out to be tachyonic for all values of the radial variable. It is shown for a sample value 
of the radial variable in Figure [H 

To quadratic order in fluctuations the action (lA.SP takes the form0 



Jmg V-^ 



i(ie-''(l-a=)-l)(«' + /3')). (A.9) 
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A prime denotes a derivative with respect to the radial variable. 
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Figure 9: This plot shows the potential for fluctuations around the cylinder solution in the 
CVMN background for a generic value of the radial variable. 



Clearly this quadratic action has an 50(2) symmetry acting on the fluctuations a and /3, 
so that only 6{t, p) = a/ «(t, py + (3{t, pY has a potential and we can focus our attention on 
this mode. Its equation of motion takes the form 

6 = 6" + 2<p'6' + (1 - \e~^^{l - a^))6 . (A.IO) 

As usual, this can be solved by separation of variables. Introducing 6(t,p) = T{t)e~'^^{p) 
we find the following system of ordinary differential equations 

f = -M^T, (A.ll) 
+ + + \e'^^{\ - a2)) ^ = . (A.12) 

We see that there will be a tachyon, i.e. an exponentially growing mode, if the second 
equation has normalizable solutions with negative eigenvalues!^ This is now just a one- 
dimensional quantum mechanics exercise with a potential 

which is shown in Figure [TOl 



^^The condition of normalizability is imposed on us from the requirement that the energy in a fluctuation 
should be finite. 
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Figure 10: This plot shows the potential (1A.13P of the Schrodinger problem (lA.12p . 



Clearly, since the potential is positive definite, there are no eigenfunctions with negative 
eigenvalues and hence no tachyons implying that the cylinder solution is indeed stable in 
the CVMN background. In fact, using the usual norm we are used to from quantum me- 
chanics0 there are no normalizable eigenfunctions at all, implying a continuous spectrum of 
fiuctuations. This might be somewhat surprising since one might expect these fiuctuations to 
correspond to part of the meson spectrum. We might wonder if one should in fact cut off the 
integral at the value of the radial variable at which the dilaton becomes too large to be able 
to trust supergravity, which would make the spectrum discrete. Fortunately, independent 
of this issue, the positivity of the potential already implies stability. Together with kappa 
symmetry this tells us that the cylinder solution is indeed supersymmetric in the CVMN 
background and it seems reasonable to smear it over the transverse directions. One might 
be worried that the smearing could introduce an additional instability but we will see in the 
next subsection that this is not the case. We attribute this to the fact that the fiber the 
branes wrap is non-trivially fibered over the base space. 

We should point out that we have chosen the orientation of the brane so that it is kappa 
symmetric as well. If one reverses the orientation of the brane the last term in the potential 
changes sign and the potential is in fact no longer positive definite and does admit a bound 
state with negative energy corresponding to a tachyon. 



A. 2. Stability of the Cylinder Solution in the Backreacted Background 

Now that we have convinced ourselves that a single probe brane in the CVMN background 
embedded as described above is supersymmetric it makes sense to try to smear it and take the 
backreaction into account. To make sure that the smearing does not spoil the supersymmetry 
we check that an additional probe D5-brane is still stable. In Appendix [B] we check that it 

^^To analyze stability in gravitational systems with a timelike Killing vector one checks that the energy 
is positive for all fluctuations for which it converges. The usual quantum mechanics norm in our case is 
equivalent to the norm induced by the energy functional. 
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is also kappa symmetric. Again we will see that the only possibly dangerous modes are the 
modes corresponding to the brane slipping off the three-sphere. The worldvolume action is 
of course still given by (1A.5I) . The six- form in our case is given by 

C(6) = ^e^'^dx^ A dx^ A dx^ A dx^ A dp A {u^ - A^) . (A. 14) 

For our embedding we find that this leads to a potential for the fluctuations given by 



V{a, l3) = e^'^ ye^ cos a cos /Jy e^3(sin^ a cos^ /? + sin^ /?) + e^'^ cos^ a cos^ {3 

- cos^ a cos^ 13^ . (A. 15) 

Again this depends on the radial coordinate via the dilaton as well as through the functions 
k and g. For all values of the radial variable for which the supergravity approximation can 
be trusted this is no longer tachyonic. The presence of the flavor branes has stabilized this 
mode@ It is shown for a sample value of the radial variable for the solution in the decoupling 
limit with = 1 in Figure [TTl 




Figure 11: This plot shows the potential for fluctuations around the cylinder solution in the 
backreacted background for a generic value of the radial variable. 

From this we can already guess that the configuration will be stable but we will be careful 

^^This statement clearly depends on the orientation of the probe brane. We have checked that for the 
kappa symmetric configuration it is stabilized while for the opposite orientation it is tachyonic. 
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and complete our analysis. To quadratic order in fluctuations the action (lA.SP takes the 
fornix 

-i(e2^-2(e^^-Ar,))(a^ + /32)^. (A.16) 

Again, this action has an 5*0(2) symmetry acting on the fluctuations a and /5, so that only 
5{t, p) = \/o;{t, pY + (3{t, pY has a potential and again we focus our attention on this mode. 
Its equation of motion takes the form 

8 = - e-2^-2^'(e2^ - 2(6^'^ - N,))5 . (A.17) 

Introducing 5{t, p(x)) = T{t)e~'^~^~^^'^'^{x) and defining a new radial coordinate x via dx = 
e^dp we find the following system of ordinary differential equations 

f = -M^T, (A.18) 
+ ((I + g' + (p'y + (f + g" + (p") + e^2fe-29(e29 - 2{e^'' - N,))) ^ = M^^A.19) 

where the prime now indicates a derivative with respect to x. As before a tachyon corresponds 
to a normalizable solution of the second equation with negative eigenvalue. Again we have 
reduced the problem to a one- dimensional quantum mechanics exercise with a potential 

V{x) = (I + ^' + '^') ' + (y + 9" + 0") + e-''-'^{e'^ - 2{e'' - iV,)) , (A.20) 
which is shown in Figure [T2l 

Somewhat surprisingly, the potential is no longer positive definite and it even seems that it 
might admit bound states with negative energy eigenvalues. However, we have to remember 
that our supergravity solution is singular and can only be trusted far enough from the 
singularity. So we argue that one should impose a hard cutoff, i.e. introduce an infinitely 
high wall, at the value of the radial variable below which supergravity breaks down. The 
value of the radial variable below which supergravity breaks down for this particular case 
is indicated by the dashed line. With a cutoff the potential will once again not lead to any 
eigenfunctions with negative eigenvalues and we conclude that tachyons are absent even once 
we have smeared the branes and taken the backreaction into account. The same remarks 
about the norm that we made for the CVMN background apply for our solutions as well 
but this issue once again does not alter the conclusion that the configuration is stable under 
small fluctuations. 

^^A prime denotes a derivative with respect to the radial variable, p. 
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Figure 12: This plot shows the potential flA.201) of the Schrodinger problem flA.lQp for the 
near horizon geometry with po = 1. 

B. Appendix: Derivation of BPS Equations and Kappa 
Symmetry 

In this appendix we derive the BPS equations of type JIB supergravity for our ansatz fol- 
lowing [T3],p^@ We then derive the kappa symmetry condition [lT],[l2j for a single flavor 
probe brane localized in the 6, ip, 6, (p directions. Since our background is a solution to the 
BPS equations, the stability of the branes shown in Appendix El together with the kappa 
symmetry implies that the system is supersymmetric |28j . 

B.l. BPS Equations 

Recall that the RR 3-form field strength in the vielbein basis (I2.10p is 



+ N^e-^^-^-^-^a{e^ A + A e^) A . (B.l) 



^°We are considering the special case h{p) = 0, a; = 2 of Ji)J and use the rescaled radial coordinate 
dp = e~^dr . 
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The non-vanishing components of the spin connections are 



LU 



UJ 



-00^^ = -a! e-f eP 
4 



00 



23 



4 

+ ae~^~^ cosh(fc - 5f)e'' 



^13 _ g-/+fc-2<;,2 



e^ + ae-^-''cosh(A;-c/)e'^ 
1 



Co" 



3ip 



lip 



,26* 



asinh(/i; - 5()e ^e' , 



2p 



'e-f-h+g-k LP 
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^_g-/+fc-2<; ^ 2e--^-'=)e3 - e-^-^ cot ^e"^ , 



l)e^ + ae--^-'^sinh(A;- c/)e^ 
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__i^^-2h-f+k^^2 _ ^ gij^j^^^ _ ^y2 _ 



The supersymmetry variations for the dilatino A and the gravitino tl)^ in the Einstein frame 
are given by |13] 



1 



62 



4 yb 



(B.2) 



where A, and e are complex Weyl spinors. 

If we write e = e+ + e_ such that ie* = e+ — e_ and then solve the variations for e+, it is 
easy to verify that this implies that they cannot be solved for e_ . So we can drop the index 
and impose that ie* = e. 

Furthermore, e = e(p) due to the symmetries of our ansatz for the metric. 
The spinor e also has to satisfy the condition that 



rfl„e — Fioe. 



fB.3) 



This follows for example from demanding that the Q dependent part of 5'?/'<^ vanishes and 
corresponds to demanding that the color branes wrap a supersymmetric 2-cycle. 
The dilatino variation is 



9(^3 



+ N,e-''f-''-^~''a{Te2 + T^i)T,]e. 



Comparing this with the components of the gravitino variation we find that 
Using this and multiplying the equation above with 2e-^ri23/0' we have 



(B.4) 
= 4/. 



(B.5) 
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where 

^=^(e-^-ie-V-l)), (B.6) 
N 

B = ^-ae-^~^. (B.7) 

(f)' 

From the conditions (rpi23)^e = e and {rpi23, r<^2} = it follows that 

A'' + B^ = 1. (B.8) 
Therefore, we can write A = cos a, B = sin a. The projection condition (IB.Sp then becomes 

rpi23e = e^^-^e. (B.9) 
In terms of eo defined by e = e~^^'^^eo this projection conditions is simply 

rpi23eo = eo- (B.IO) 
From the p component of the gravitino variation we find 

d,e-^-]a'e-^+^T^2e = 0, (B.ll) 
o 4 

which gives 

dpeo - jT^.eo - ^ - \a'e-''^^T^,eo = 0. (B.12) 
Solving the terms with and without r<^2 independently we have 

a' = -^a'e-^^' and eo = e'^/^rj, (B.13) 

where ?7 is a constant spinor satisfying the same projections as cq. 
So to summarize, our Killing spinor e is of the form 

e = e'^^'^'^esT] , (B.14) 

where ?7 is a constant spinor satisfying the projections 

^pi23V = T], Te^T] = TuV, if]* = V 

and a is obtained by integrating 

ry' = - . 
2 



From the remaining variations we obtain two equations each, one from the part proportional 
to e and the other from terms proportional to r^2e- 
From Sipxt^ = or = we get 

/' = ^Ae-^^-—Ae-^\a^-l)--aBe-''~^N,, (B.15) 
4 16 4 

= Aa + Be^-' -^Bef^-^a^ -1). (B.16) 
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From Sipe = or equivalently from 5'^l)^p = and using the previous expressions we get 
h' = -Bae-''+'' (^cosh(fc - g) ~ y""'"'') ~ \-Ae-^''^^\a^ - 1) 

+^Ae-^\a^ - 1) , (B.17) 
a' = -AAae-3+^ cosh(A: - g) + Be-3-^+'^^{a^ - 1) - 2NJ3e^-^^ 

-^Be~''~'{a^ - 1). (B.18) 

From Sipi = or 5il)2 = we get 

g' = Ae'^-^^ -Bae-^+Hmh{k-g) + ^Bae-^-3 -N^Ae-^\ (B.19) 

= Aa + Be^-3 -^Be^-^{a^ -I). (B.20) 

Note that (lR20ll is the same as (fRTHll . 
From 5ipz = and (1B.16P one finds 

k' = -Ae^''-^^ + 2A + ^Ae-^''+^\a^ -l) + 2aBe-''+''smh{k-g) 

-N,Ae~^' + ^Ae-^\a^ - 1) + N.Bae-'"^ , (B.21) 

= 2B + aAe3-\ (B.22) 
It is straightforward to check that for a(p) = ^^^2p these equations imply 

e''-9 = - tanh 2p, A = tanh 2p, B = 1 — , (B.23) 

2 cosh 2p 

and 

^^e'^+s = e2fc _ ^3 24) 

dpk = -(e^'^ + N,)e-''-3 + 2 coth 2p, (B.25) 

dpf = ^e-^-^ (B.26) 

B.2. Kappa Symmetry 

Now we follow [13j and show that a single probe fiavor brane in the backreacted background 
is kappa symmetric. We take the fiavor brane to extend along the x^, p and directions and 
at fixed 9, ip, 9, (p. For e = ie* and in the absence of a worldvolume gauge field and 5-field 
we have to show that 

r.e = e, (B.27) 

where 

1 1 rn\ m% 

= 77 / 7ml. ..m6- (B.28) 
6! y/-5({6) 
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g-mn = dmX^dnX^ g^j_y and 7^ = dmX^e'^^a are the puUback of the metric to the worldvolume 
of the brane and the induced Dirac matrices on the worldvolume of the brane, respectively. 
For our embedding we can choose x^, and p as worldvolume coordinates. This gives 

7.. = e^r.., 7^ = ^13, 7p = e^+%, = -e'f+^'. (B.29) 

So we need to show that 

r^e = r2.0a.ij.2^3 3pe = e. (B.30) 

Since e is a spinor of definite chirality of type IIB supergravity we know that it satisfies 
^x0x^x^x3pe>fi23(^ = e- From multiplying ([RS]) by T^o^i^2^3e^ we get 

T^o^i^2^3{cosaT0^ + sina;r5i2)e = e, (B.31) 

and from multiplying fIB.Sp with Tq^ one can show that 

(cos aTg^ + sin 01512)6 = T^pe. (B.32) 

So we find that 

r^e = ra.0a,i^2^33pe = e , (B.33) 

which means that the probe brane is kappa symmetric. Note that this derivation does not 
depend on the angular coordinates and is therefore true independent of the brane's position. 



C. Appendix: Analytic Solutions 

In this appendix we derive the analytic solutions given above in Section 13.21 compare them 
to the numerical solutions and discuss some of their properties. 



C.l. The Solution in Region I 

This region is characterized by e"^^ <^ Nc- One might simply drop e"^^ from the equations 
(13.41) but we will take a more systematic approach and treat e^'^ as a small perturbation. If 
we define 

e^+^ = iV,(7 + 57), (C.l) 
e^'' = NM, (C.2) 

and keep only terms up to first order in small quantities, the equations become 



dpi = -1, 
(9p(57 = 5k, 

2 

dp\n6k = -- + 4coth2p. 

7 



(C.3) 
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This set of equations has the the solution 

e^+VA^c =Po - P + ^^^2 (i (1 + 8(P - P^f sinh 4p) (C.4a) 

g 

- (P - Po) cosh4p - po + -p(p^ - 3ppo + 3po; 

e-/iV^ = f!i^fl-iL^V. (C.4b) 



2p/ V Po 
To zeroth order in e^^ the result is 



/N, = po-p, (C.5a) 
e-/iV.= f^fl-i^')V. (C.5b) 



V 2p/ V Po 

Solutions Contained entirely in Region I 

In order for solutions to be contained entirely in Region I we have to impose that e"^^ <^ 
for all values of p. This is satisfied if for all p 



2p7 V Po 

For these solutions e^'^ has a maximum and it is enough to impose that the value at the 
maximum be much less than Nc- The maximum is located at a value p given by 



' ^ 

dp 



smh '^P f^_p\\ _1 ^2(po - p) cosh 2p - sinh 2p) . (C.7) 



So p is given by the solution of 



2p/ V Po// 2poP7 



2(po - p) = tanh2p. (C. 



From this it is clear, that p will only depend on po and not on pj and we can rewrite the 
equation (]C.6P as a relation between pi and po 

1 . - ~ , ^ p(po) \ 



p/ > -sinh2p(po) 1 = pi,min{po)- (C.9) 

Po / 



This can be solved numerically. The result is shown in Figure [131 
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Figure 13: The shaded region in this plot shows the region in po ~ Pi parameter space for 
which the solutions are contained entirely in Region I. The solid line shows the function 
Pl,min (Po)- 



Comparison of Numerical and Analytic Solution in Region I 



To get an idea of how good the analytic solutions given in (]C.4p are, we compare them to 



the exact numerical solutions for = 1 m Figure [T3] for three different values of p/ given 
by 2, 4, and 8 times the bound pi,min according to ( 1C.9I1 . We see that the further we move 
inside the shaded region of Figure [13] the better the approximation becomes. We also see 
that for values of the parameters po and p/ near the border we should keep the next order 
in our expansion for e^^ and that away from the border even the approximation fIC.SP is in 
excellent agreement with the numerical results. 

C.2. The Solution in Region II 

This region is characterized by e^^^ ^ N^.-, e^^ . Under these assumptions the equations (13.41) 
become 

' (CIO) 
Opk = 2 coth2p. 

These equations have the solution 



e-/^. . ,c.n, 

.-/iV. = (C-12) 
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Figure 14: Comparison of the exact numerical solution shown as a solid line and the approx- 
imate solution according to flC.4l) shown as a dashed line for pq = 1 and pj = 2pi jnin(Po)) 
Pi = 4pi jnm(Po)) and pj = 8pi^min(po)) in the top, middle, and bottom row, respectively. 

Clearly these solutions will only be valid as long as 

p sinh4p ^ /sinh2pV /^-,o\ 
8p/^ 32pi^ \ 2pi J 
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C.3. The Solution in Region III 

This region is characterized by e^^ ^ Nc. The equations in this region become 

p4fc (C.14) 
a^e'^ = -2-— + 46^^= coth 2p. 



These have the solution 



eS+^Nc = A{8B- 12p + 3 sinh4p) ^ , (C.15a) 

e^ViVc = ^A^I^ ^. (c.lSb) 

(85- 12p + 3sinh4p)3 

Here A and B are so far arbitrary integration constants. We will see how to relate them to 
more useful quantities below. 

Solutions Contained entirely in Region III 

For solutions to be contained entirely in Region III the condition e^'^ ^ Nf. has to hold for 
all p, in particular for p ^ 0. This implies that solutions contained entirely in Region III 
must have 

B = 0. (C.16) 



So the solutions flC.lSP become 

e^+'^/N^ = 3^A(sinh4p-4p)^ (C.17) 

e-/iV. = J^^'^^'^P ^ . (C.18) 
33 (sinh 4p — 4p)3 



After identifying 



4 
£3 



A=-r-T, (C.19) 

2333 



these can be seen to correspond to the deformed conifold [29j 



.^^VA^c = ^^(^^"h4f-4^)\ (C.20a) 

23 



— 1 2 

2k /jKT 2 4 23 sinh 2p 



e-/iV, = -e^— — f-j. (C.20b) 

•J (smh4p — 4pj3 



In order for this to be contained entirely in Region III we require 

e»Qy^l. (C.21) 

Comparison of Numerical and Analytic Solution in Region III 

A comparison of these approximate solutions to the exact numerical solutions in this region 
are shown in Figure [151 As was to be expected the approximate and numerical solution 
agree better and better the larger the deformation parameter e. 
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Figure 15: Comparison of the exact numerical solution shown as a solid line and the approx- 
imate solution according to (1C.20P shown as a dashed line for e = 10, e = 100, and e = 1000, 
in the top, middle, and bottom row, respectively. 



C.4. Matching Solutions in Regions II and III 

The matching of the solutions in Regions II and III cannot be done in general but we can 
do so in two limiting cases, namely p/ -C 1 and 3> 1. 

Matching for pj <^ 1 

In the case of small pi e^'^^ approaches its minimum value for small p and for small enough 
Pi the matching can be done for values of p that are small compared to unity but large 
compared to pj. For p the solutions in Region II and III take the following form: 
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Region II 



o3 



e^+V^e ^ Po + ^, (C.22) 

e'VTVc ^ (C.23) 

In particular this implies that we need po 1 to ensure the validity of the approximations 
we used to find this solution. 



Region III 



This implies 



2 \2pi 

So the solution can be written as 



.g+hj^^ _^ (C.24) 
e''/N, - ^p'. (C.25) 



2 

A=- [P^] ' and B = 4poPi. (C.26) 



and similarly 



e^V^^c = <! , , 2 . ^ ^ , , (C.28) 



sinli 2p 
2pi 



where is the value of p for which we match the two solutions. For our approximations 
to be valid we require po ^ 1 as well as p/ ^ p* ^ ^/PoPi- Wc can choose a p* such that 

p* <^ (PoP/)^- Introducing a small additional error especially for small p for p/ <C 1 we can 
write the solution as follows for all values of p 

e-V/V. = p„ (l - ^ (3p - 5 sinh 4p) ) ' . (C.29a) 



/ sinh2p\ / 1/3 



\ 2pi J \ Spopi 



e-/iV,= (^^) (l-^TT^(3p--sinh4p) ) . (C.29b) 



Matching for p/ » 1 



For Pi ^ 1 the function e^^^ takes its minimum for p 3> 1 and so the matching will also 
be done for some value p* » 1. In this limit the solutions in Regions II and III take the form: 
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Region II 



e'+ViVc - Po + ^, (C.30) 

(C.31) 

In particular this implies that we need po 3> 3> 1 to ensure the validity of the approxi- 
mations we used to find this solution. 



Region III 



e^+ViVc ^ A(8B+^e^py , (C.32) 

e^VA^c - ^ (C.33) 
(85 + fe^'') 3 



Matching these then leaves us with the same conditions as before 



2 

A=^ (^) ' and 5 = 4poP/. (C.34) 



So again the solution can be written as 



1 p I sinh4p . 



^0 (} ~ 8^ (3P - I sinh 4p) j p > 



and similarly 



e^'/N^={ . ^ . , (C.36) 

(^) (l-8^(3p-!-nh4p)) ' p>p. 

where now the value of has to be chosen so that it satisfies 

^ ln4pj < < ^ ln4^pj. (C.37) 

Clearly this is only possible if we take po ^ 1 but this is consistent with what we have 
assumed above. If we take po ^ p* this simplifies further and we can, again introducing a 
small additional error, write the solution in closed form 



=»^ViVe = P„(l + ^V-A (C.38a) 



e^ViVe = ^ 1 + ; o • (C.38b) 



2 



/sinh2pV/ 2>e^P ^ ^ 
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Figure 16: Comparison of the exact numerical solution shown as a solid line and the ap- 
proximate solution according to ( 1C.29P shown as a dashed line for po = 8 and pi = 10"^, 
Pi = 10~^, and pi = 10~^, in the top, middle, and bottom row, respectively. 

Comparison of Numerical and Analytic Solutions 



The approximation (]C.29|) is compared to the numerical solution for po = 8? for pi = 10 



Pj = 10 ^, and p/ = 10 ^ in Figure [T6l As can be seen the approximation works reasonably 
well and the accuracy increases as p/ decreases. The approximation (]C.38P is compared to 



the numerical solution for po = 20, and for p/ = 5, p/ = 10, and p/ = 20 in Figure [T71 Again 
the approximation works reasonably well and the accuracy increases as p/ increases. 
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Figure 17: Comparison of the exact numerical solution shown as a solid line and the approx- 
imate solution according to flC.38|) shown as a dashed line for po = 20 and pi = 5, p/ = 10, 
and pi = 20, in the top, middle, and bottom row, respectively. 
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